:arXiv:1507.03777v2 [cond-mat.mtrl-sci] 4 Dec 2015 


Nitrogen as the best interstitial dopant among X=B, C, N, O and F 
for strong permanent magnet NdFenTiX: First-principles study 


Yosuke Harashima,^’^ Kiyoyuki Terakura,^’"'^ Hiori Kino,^’^ Shoji Ishibashi/ and Takashi Miyake^’^ 
^Nanomaterials Research Institute, AIST, Tsukuba, Ibaraki 305-8568, Japan 
‘^MANA, National Institute for Materials Science, Tsukuba, Ibaraki 305-0044, Japan 
^ESICMM, National Institute for Materials Science, Tsukuba, Ibaraki 305-0047, Japan 
^National Institute for Materials Science, Tsukuba, Ibaraki 305-0047, Japan 
(Dated: December 7, 2015) 

We study magnetic properties of NdFenTiX, where X=B, C, N, O, and F, by using first-principles 
calculations based on density functional theory. Its parent compound NdFenTi has the ThMni 2 
structure, which has the symmetry of space group li/mmm, No. 139. The magnetization increases 
by doping B, C, N, O, and F at the 2b site of the ThMni 2 structure. The amount of the increase 
is larger for X=N, O, F than for X=B, C. On the other hand, the crystal field parameter (r^)A^, 
which controls the axial magnetic anisotropy of the Nd 4/ magnetic moment, depends differently 
on the dopant. With increase of the atomic number from X=B, (r^)A 2 increases, takes a maximum 
value for X=N, and then turns to decrease. This suggests that in NdFenTiX, nitrogen is the most 
appropriate dopant among B, C, N, O, and F for permanent magnets in terms of magnetization 
and anisotropy. The above calculated properties are explained based on the detailed analysis of the 
electronic structures of NdFenTiX. 


I. INTRODUCTION 

Large remanent magnetization and high coercivity are 
the two fundamental requirements for high performance 
permanent magnets. The remanent magnetization is de¬ 
termined by spontaneous spin and orbital moments of 
a material. The coercivity is strongly correlated to the 
magnetocrystalline anisotropy energy. Extensive studies 
have been devoted to the search for higher-performance 
permanent magnets. The strongest permanent magnet 
to date is based on Nd 2 Fei 4 B, which has magnetization 
of 1.85 T at 4.2 K, magnetocrystalline anisotropy field of 
67 kOe at room temperature, and Curie temperature of 
586 K.- NdFeiiTiN (Fig.[T|) is also a strong magnet com¬ 
pound, although its maximum energy product is smaller 
than Nd 2 Fei 4 B. The interstitial-nitrogenated NdFenTi 
has magnetization of 1.476 T at 1.5 K, magnetocrys¬ 
talline anisotropy field of 80 kOe at room temperature, 
and Curie temperature of 740 
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FIG. 1. (Color online) Crystal structure of NdFenTiN. Ti 
occupies one of the 8i sites. 


NdFenTi has the ThMni 2 structure (space group 
lA/mmm, No. 139). Ideally all of the Mn sites are occu¬ 
pied by Fe. However, NdFei 2 is thermodynamically un¬ 
stable, and some of Fe have to be substituted by another 
element to stabilize the bulk phase (NdFei 2 -xT/a;). For 
example, Ti, V, Cr, Mn, Mo, W, Al, and Si are known to 
serve as such substitutional elements However, as 

these substitutional elements reduce the magnetization^, 
the concentration of the substitutional atoms should be 
as small as possible. The concentration range of the sub¬ 
stitution, X, depends on the element. In this regard, ti¬ 
tanium is a favorable element because it stabilizes the 
material with small x (« 1). 

To optimize the performance of a strong magnet, 
we can utilize even another degree of freedom, i.e., 
the interstitial dopant. In fact, magnetic properties of 
NdFei 2 -xAIa, are significantly controlled by interstitial 
doping of a typical element. It is experimentally observed 
that the interstitial nitrogenation enhances the mag¬ 
netic moment from 21.273 (NdFenTi) to 23.218 /is/f.u. 
(NdFeiiTiNo.s);^ The Curie temperature (Tc) also rises 
from 570 K in NdFenTi to 740 K in its nitridej^i^ Fur¬ 
thermore, strong uniaxial magnetocrystalline anisotropy 
is induced by the nitrogenationj^ii Another possible 
interstitial dopant is carbon<^“— Carbon doping at an 
interstitial site leads to increase of Tq and uniaxial 
anisotropy in NdFei 2 -xAl 2 ; compounds. However, nitro¬ 
genation is more preferable than carbonation in terms of 
magnetization and Curie temperature. In order to de¬ 
sign better strong magnets, the cause of the difference 
between N and C has to be clarified and the possibility 
of better dopants has to be searched for. 

In previous papersjiid^ we have theoretically studied 
the effects of Ti substitution and nitrogenation by com¬ 
paring NdFei 2 , NdFenTi and NdFenTiN using exper¬ 
imental structures. We found that Ti substitution in 
NdFei 2 reduces the magnetic moment more than sub- 
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tracting the local magnetic moment at the substituted 
Fe site, and enhances the uniaxial anisotropy slightly. It 
was also shown that interstitial nitrogenation increases 
the magnetic moment by 2.75 /rs/f.u. in NdFenTi. 
The magnetic anisotropy of NdFenTi and its nitride was 
studied by using the crystal field parameter We 

found that the interstitial nitrogenation substantially in¬ 
creases ir^)A 2 . This implies that the uniaxial anisotropy 
is enhanced by the nitrogenation in NdFenTi. Therefore, 
the interstitial nitrogenation works preferably in terms of 
magnetization and magnetocrystalline anisotropy. 

Stimulated by our work, the experimental group tried 
to synthesize a new strong magnet compound NdFei 2 N 
and succeeded in forming a film recently.i^^ It has been 
shown that the film sample has better intrinsic magnetic 
properties than Nd 2 Fei 4 B, although a bulk sample is 
still difficult to synthesize. This experimental work, in 
turn, motivated us to revisit NdFenTiN and related com¬ 
pounds in order to search for stronger permanent magnet 
compounds. 

In the present paper, we perform detailed electronic 
structure calculations for NdFenTiX with X the inter¬ 
stitial impurity of typical elements B, C, N, O, and F. We 
analyze in detail the electronic structures, paying partic¬ 
ular attention to the X-2p states. Then, based on the 
analysis we explain the variations in structural proper¬ 
ties, magnetization, and magnetocrystalline anisotropy 
caused by the interstitial dopant X. 


II. CALCULATION METHODS 

The calculations are carried out by using the first 
principles code QMAS (Quantum MAterials Simulator)ii 
which is based on density functional theor} ^^^^ and the 
projector augmented-wave (PAW) methodJ^d^ For the 
exchange-correlation energy functional, the generalized 
gradient approximation (GGA) is usedj^^ The 8x8x8 
k points are sampled and the cutoff energy for the plane 
wave basis is set to 40.0 Ry. The lattice constants and 
inner coordinates are optimized. NdFenTiA has 4/ or¬ 
bitals that are strongly localized at the Nd sites. In the 
present study, the /-states of Nd are treated as open-core 
states, in which the hybridization with other orbitals is 
neglected completely and the atomic physics is applied to 
the Nd-4/ states. Following Hund’s first rule, we assume 
that Nd has three /-electrons with the full spin polariza¬ 
tion. The self-consistent calculation makes the 4/ spins 
antiparallel to the Fe spins. 

The orbital magnetic moment is not treated in our self- 
consistent calculations, but that of the Nd-4/ electrons 
is included in the total magnetic moment as follows. We 
assume that the Nd-4/ electrons yield the local mag¬ 
netic moment of gjJ, where gj is the Lande g-factor, 
and J = 9/2 is the total angular momentum for the Nd- 
4/ electrons. (For Nd-4/ electrons, J = \L — S'], the 
orbital angular momentum L = 6, and the spin angular 
momentum S = 3/2 are given by Hund’s rules.) The to¬ 


tal magnetic moment of the whole material is estimated 
by adding gjJ=3.273 pB to the spin magnetic moment. 
This spin magnetic moment is the contribution from the 
valence electrons other than the Nd-4/ electrons. Fe has 
only small spin-orbit interaction compared with Nd; thus 
the spin-orbit interaction on Fe is neglected. 

The magnetocrystalline anisotropy energy is expressed 
in terms of a first-order coefficient Ki as 

E{e) K. Ki sin^ 9 , (1) 


where 9 is the rotation angle of the magnetization mea¬ 
sured from the c axis. We focus on the contribution of the 
rare-earth 4/ electrons. It can be evaluated by the inter¬ 
action energy between the 4/ electrons and the surround¬ 
ing crystalline electric field. Because of Hund’s rules, the 
electron distribution of the Nd-4/ states deviates sig¬ 
nificantly from the spherical symmetry. The shape of 
the electron density distribution is uniquely related to J 
which is parallel to the 4/ magnetic moment. The non- 
spherical electron density distribution couples with the 
crystal-electric field leading to the magnetic anisotropy. 
Following the crystal-field theory, Ki can be expressed 
by using the crystal field parameter A 2 as 

Ai = -3J(J-i)aj(r2)A°nfl. (2) 

Here, aj is the Stevens factor of the rare-earth atoms, 
which is determined for each rare-earth element. For 
Nd^’*', aj = —7/1089. ub is the Nd concentration. We 
evaluate {r^)A^ by the equation 


(r-')A)" = Fr [ " Wr (r) (r) dr . 

Jo 
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F[^ is a prefactor of the real spherical harmonics ZF and 
its explicit expressions can be found in, e.g.. Ref. 1^. IF™ 
is the effective potential at the Nd site expanded by ZJ^. 
(j) is the radial function of the Nd-4/ orbital, which is 
obtained in GGA with the self interaction correction. In 
Eq. (21) it is implicitly expected that (j) is well localized, 
and we neglect the contribution from a tail of (f> in Eq. m 
by introducing a cutoff radius Vc- The cutoff radius is de¬ 
termined so that the volume of the sphere is equal to that 
of the Bader region i^^i^^ (See Fig. SI of the Supplemen¬ 
tal Material^^ for the actual value of the atomic sphere 
for each X. The atomic radius for NdFenTiA is also 
used for NdFenTiEjf, which denotes NdFenTi given by 
removing X from NdFenTiA with other atoms hxed at 
their positions in the structure optimized.) 


III. RESULTS 

From the total energy of NdFenTi, we found that Ti 
substitution is more stable at the 8i site than the 8/ 
and 8/ sites by 0.78 and 0.51 eV/f.u., respectively, which 
agrees with experimental indicationi^i^ (Even though 
the symmetry of the ThMni 2 structure is broken by the 
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TABLE 1. Optimized lattice constants of NdFenTiX in nnits 
of A. Due to Ti substitution, the system is no longer tetrag¬ 
onal but orthorhombic {a ^ b). 



a 

b 

c 

volume 

X^B 

8.521 

8.593 

4.943 

180.9 

C 

8.495 

8.572 

4.914 

178.9 

N 

8.537 

8.618 

4.880 

179.5 

0 

8.658 

8.704 

4.786 

180.3 

F 

8.812 

8.830 

4.745 

184.6 

empty 

8.553 

8.568 

4.701 

172.3 



1.78 

1.76 

1.74 

1.72 

1.70 

1.68 

1.66 

1.64 

1.62 

1.60 

1.58 

1.56 
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FIG. 3. (Color online) X dependence of {r'^)A% and K\. The 
filled red circles and empty black circles correspond to 
of NdFeiiTiX and NdFenTiEx, respectively. The latter is 
the system without X keeping the other atoms fixed at the 
positions of the former. The lattice constant and inner co¬ 
ordinate are optimized for NdFenTiX. (The parameters are 
shown in Tables [J and M) The structural effect shown by 
empty black circles is not dominant in the dopant dependence 
of (r^)j42. The details are discussed in Sec. lIVDl In addition, 
Ki is estimated from {r^)A!^ by using Eq. and shown as 
the blue squares. The red circles and blue squares do not scale 
exactly because of the volume effect through n_R in Eq. (O- 


FIG. 2. (Color online) Total magnetic moment and magneti¬ 
zation of NdFeiiTiX (A=B, C, N, O, F) in units of /rs/f.u. 
(left axis) and tesla (right axis). Theoretically optimized 
lattice constants and inner coordinates are used. In all the 
dopants, the magnetization is larger than that of NdFenTi, 
1.631 T, shown as a blue dotted line (corresponding magnetic 
moment is 24.10 /rs/f.u. shown as a red broken line). 


Ti substitution, we use the notation of ThMni 2 for the 
structure of NdFenTiAT.) Hereafter, Ti is put at the 81 
site and X is inserted at the midpoint between Nd atoms 
along the c axis as shown in Fig. [T] The optimized lattice 
constants and fractional coordinates for NdFenTi and 
NdFeiiTiAi are shown in Tables U and m respectively, 
where ’’empty” denotes NdFenTi. The interstitial X ex¬ 
pands the volume by 5% for X=B, 4% for C, 4% for N, 
5% for O, and 7% for F. 

In the present paper, we distinguish between ’’mag¬ 
netic moment” m and ’’magnetization” fj,oM. ’’mag¬ 
netic moment” is expressed in units of the Bohr mag¬ 
neton fj,B per atom or per formula unit, while ’’magne¬ 
tization” is in units of Tesla which is estimated from m, 
volume, and the vacuum permeability /tq. Ki is evalu¬ 
ated from the obtained by using Eq. d^]). We also 

estimate the magnetocrystalline anisotropy field fioHa 
(= 2^oKi/ fioM). The calculated magnetic properties 
are summarized in Table mi 

The total magnetic moments are calculated to be 24.10 


/is/f.u. for NdFeiiTi and 26.84 /is/f.u. for NdFenTiN 
with the structure optimization for each system. They 
correspond to the magnetization of 1.631 T and 1.743 
T, respectively. In the previous paper, we reported 
the spin magnetic moment for the experimental lattice 
parametersiii The results were 19.99 and 20.97 /rs/f.u. 
for NdFeiiTi and NdFenTiN, respectively. Here, the 
spin magnetic moments of the valence electrons other 
than that of the Nd-4/ electrons are 22.99 and 23.97 
/is/f.u. The corresponding total magnetic moments can 
be estimated by adding gjJ=3.273 fis to these spin mag¬ 
netic moment as 26.26 ^s/f.u. (1.697 T) and 27.24 

/rs/f.u. (1.731 T), respectively. The nitrogenation en¬ 
hances the magnetization for both optimized and experi¬ 
mental lattice parameters, and the increment is larger in 
the former case (0.112 T) than the latter (0.034 T). Ex¬ 
perimentally, the magnetic moment is 21.273 /rs/f.u. in 
NdFeiiTi and 23.218 /rs/f.u. in NdFenTiNo.sr which 
correspond to 1.375 T and 1.476 T. The increment of the 
magnetic moment and magnetization are 3.89 g,B and 
0.202 T per nitrogen, respectively. In Ref.0, on the other 
hand, the magnetizations of NdFenTi and NdFenTiNi.s 
are reported as 1.70 T and 1.84 T, respectively; thus the 
magnetization difference is 0.093 T per nitrogen. These 
experimental results are in reasonable agreement with 
our calculation for the optimized lattice parameters. 

Figure [2] shows the magnetic moment and magnetiza¬ 
tion for X=B, C, N, O, and F. The magnetic moment is 
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TABLE II. Optimized inner coordinates. The Nd atom is fixed at the origin. Even though the symmetry of the ThMni 2 
structure is broken by the Ti substitution, we use the notation Fe(8/), Fe(8i), and Fe(8j) in NdFenTiA. 



Fe(8/) 

Fe(8i) 

Fe(8i) 

A 

X=B 

( 0.257, 0.252, 0.250 ) 

( 0.257, 0.748, 0.750 ) 

( 0.757, 0.248, 0.750 ) 

( 0.757, 0.752, 0.250 ) 

( 0.374, 0.000, 0.000 ) (Ti) 

( -0.351, 0.000, 0.000 ) 

( 0.007, 0.359, 0.000 ) 

( 0.007, -0.359, 0.000 ) 

( 0.274, 0.500, 0.000 ) 

( -0.262, 0.500, 0.000 ) 

( 0.510, 0.272, 0.000 ) 

( 0.510, -0.272, 0.000 ) 

( 0.006, 0.000, 0.500 ) 

C 

( 0.256, 0.252, 0.250 ) 

( 0.256, 0.748, 0.750 ) 

( 0.756, 0.248, 0.750 ) 

( 0.756, 0.752, 0.250 ) 

( 0.374, 0.000, 0.000 ) (Ti) 

( -0.351, 0.000, 0.000 ) 

( 0.007, 0.360, 0.000 ) 

( 0.007, -0.360, 0.000 ) 

( 0.281, 0.500, 0.000 ) 

( -0.269, 0.500, 0.000 ) 

( 0.509, 0.279, 0.000 ) 

( 0.509, -0.279, 0.000 ) 

( 0.005, 0.000, 0.500 ) 

N 

( 0.256, 0.252, 0.250 ) 

( 0.256, 0.748, 0.750 ) 

( 0.756, 0.248, 0.750 ) 

( 0.756, 0.752, 0.250 ) 

( 0.373, 0.000, 0.000 ) (Ti) 

( -0.352, 0.000, 0.000 ) 

( 0.006, 0.360, 0.000 ) 

( 0.006, -0.360, 0.000 ) 

( 0.279, 0.500, 0.000 ) 

( -0.269, 0.500, 0.000 ) 

( 0.510, 0.277, 0.000 ) 

( 0.510, -0.277, 0.000 ) 

( 0.005, 0.000, 0.500 ) 

0 

( 0.256, 0.251, 0.250 ) 

( 0.256, 0.749, 0.750 ) 

( 0.756, 0.249, 0.750 ) 

( 0.756, 0.751, 0.250 ) 

( 0.373, 0.000, 0.000 ) (Ti) 

( -0.351, 0.000, 0.000 ) 

( 0.006, 0.360, 0.000 ) 

( 0.006, -0.360, 0.000 ) 

( 0.271, 0.500, 0.000 ) 

( -0.263, 0.500, 0.000 ) 

( 0.509, 0.270, 0.000 ) 

( 0.509, -0.270, 0.000 ) 

( 0.003, 0.000, 0.500 ) 

F 

( 0.257, 0.251, 0.250 ) 

( 0.257, 0.749, 0.750 ) 

( 0.757, 0.249, 0.750 ) 

( 0.757, 0.751, 0.250 ) 

( 0.373, 0.000, 0.000 ) (Ti) 

( -0.349, 0.000, 0.000 ) 

( 0.006, 0.358, 0.000 ) 

( 0.006, -0.358, 0.000 ) 

( 0.252, 0.500, 0.000 ) 

( -0.246, 0.500, 0.000 ) 

( 0.510, 0.253, 0.000 ) 

( 0.510, -0.253, 0.000 ) 

( 0.005, 0.000, 0.500 ) 

empty 

( 0.254, 0.251, 0.251 ) 

( 0.254, 0.749, 0.749 ) 

( 0.754, 0.249, 0.751 ) 

( 0.754, 0.751, 0.249 ) 

( 0.374, 0.000, 0.000 ) (Ti) 

( -0.350, 0.000, 0.000 ) 

( 0.005, 0.358, 0.000 ) 

( 0.005, -0.358, 0.000 ) 

( 0.269, 0.500, 0.000 ) 

( -0.264, 0.500, 0.000 ) 

( 0.506, 0.274, 0.000 ) 

( 0.506, -0.274, 0.000 ) 



TABLE III. Magnetic moment m [/tb/Lu.], magnetiza¬ 
tion fj,oM [T], crystal field parameter [K], magne¬ 

tocrystalline anisotropy energy coefficient Ki [MJ/m®], and 
anisotropy field fioHa [T] of NdFenTiX are shown for X =B, 
C, N, O, F, and NdFenTi. 



m 

[ mb/Lu. ] 

MoM 

[T] 

[K] 

Ki 

[ MJ/m^ ] 

[T] 

X=B 

25.51 

1.643 

-27 

-0.70 

-1.1 

C 

25.18 

1.640 

96 

2.6 

4.0 

N 

26.84 

1.743 

397 

10.6 

15.3 

0 

27.23 

1.760 

-7 

-0.2 

-0.3 

F 

27.07 

1.709 

-502 

-13.0 

-19.2 

empty 

24.10 

1.631 

-29 

-0.80 

-1.2 


increased by all the interstitial dopants. The correspond¬ 
ing magnetization is also larger than that for NdFenTi 
(1.631 T) bnt the amount of increase is partially canceled 
by the volume expansion. The increment of the magne¬ 
tization for X—B and C is small, while that for Ar=N, 
O, and F is significant with a jump between X=C and 
N. Consequently, the interstitial doping of N, O, and F 
works positively for permanent magnets than doping of 


B and C in terms of the magnetization. 

We now turn to the discussion on magnetocrystalline 
anisotropy based on the lowest crystal field parameter 
A 2 at the Nd site. Since Nd has negative Stevens factor, 
positive (negative) value of A 2 implies uniaxial (in-plane) 
anisotropy. In the case of X =N, the interstitial doping in¬ 
creases the value of {r^)A^ from —29 K to 397 K. Accord¬ 
ing to Eq. these values correspond to Ki = —0.80 
MJ/m^ and 10.6 MJ/m^, respectively. Thus, strong 
uniaxial anisotropy by nitrogenation is indicated from 
the calculation. The experimental anisotropy energy 
constant can be estimated from the magnetization and 
anisotropy field as Ki = 6.75 MJ/m^ in NdFenTiNo.s at 
1.5 K— or larger than 18 MJ/m^ in NdFeiiTiNi .5 at 4.2 
K— . If we simply assume a linear relation between dopant 
concentration and the increment of ATi, our theoretical 
estimate is 4.9 MJ/m^ for the former with N 0.5 and 16.3 
MJ/m^ for the latter with N 1 . 5 . 

The interstitial dopant dependence of is shown 

in Fig. O {r^)A 2 is —27 K for X=B. With increasing 
atomic number, {r^)A 2 increases up to A=N, then turns 
to decrease, and takes a large negative value for A=F. 
Thus, the N doping is suggested to induce the strongest 
uniaxial anisotropy among the typical elements, B, C, N, 
O,and F. 
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IV. DISCUSSION 

A. Electronic Structures 

1. Fe-3d band 



eV 


FIG. 4. (Color online) Density of states for NdFei 2 (red line) 
and bcc Fe (green line). DOS of bcc Fe has a shoulder in the 
majority-spin band at the Fermi level, and the width of DOS 
of NdFei 2 is narrower than that of bcc Fe. The origin of the 
energy is set at the Fermi level. 

As our analysis of the previous paper— and the present 
work depends on the fact that the Fermi level lies above 
the majority-spin Fe-Sd band, we first confirm this fact 
by studying the density of states (DOS). Figure |4] shows 
DOSs of both of NdFei 2 and bcc Fe. Because of the pres¬ 
ence of a shoulder of DOS for bcc Fe at the Fermi level 
in the majority-spin state, the width of the majority-spin 
Fe-3d band looks to be slightly narrower in NdFei 2 than 
that in bcc Fe. In order to make the comparison more 
quantitative, we estimated the second moment of the Fe- 
3d band DOS with the tight-binding picture using the 
neighboring Fe-Fe distances and the coordination num¬ 
bers. The inverse fifth power law was assumed for the 
interatomic distance dependence of the d-d hopping inte¬ 
gral predicted by the canonical band picture— We found 
that the second moment for NdFei 2 is 0.91 of that for bcc 
Fe being consistent with the difference in the band width 
mentioned above. This implies that the effective Fe-Fe 
distance is elongated in NdFei 2 compared with that in 
bcc Fe. The narrower d band width of NdFei 2 makes the 
majority-spin Fe-Sd band virtually filled, and is expected 
to give larger magnetic moment. The spin moment per Fe 
is 2.21 /tb in NdFei 2 and 2.18 hb in bcc Fe. This is con¬ 
sistent with the description above. However, the ratio of 
the spin moments 2.18/2.21 = 0.99 is close to 1, whereas 
that of the second moment is 0.91. In NdFei 2 there is 
also the hybridization between Nd-Sd and Fe-3d. As Nd- 
5d states are located above Fe-3d states, this additional 


hybridization pushes down the minority-spin state more 
than the majority-spin state of Fe-3d and reduces the 
magnetic moment in NdFei 2 , which explains such small 
difference in the spin moment between NdFei 2 and bcc 
Fe. The situation is not basically changed in NdFenTi 
except the appearance of the split-off Ti-3d states.— 

2. X-2s and 2p states 



FIG. 5. (Golor online) Atomic eigen levels of X-2p and -2s 
states are shown for X =B, G, N, O, and F. Fe-3d and -4s 
levels are also shown. The vacuum level is taken as the origin 
of the vertical axis. 
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FIG. 6. (Color online) The local coordinate defined with re¬ 
spect to the pair between a particular X (blue circle) and a 
particular Fe(8d) (red circle). 

In order to understand the effects of the interstitial el¬ 
ement X on magnetic properties, we have to understand 
the basic properties of the electronic structure associ¬ 
ated with X, particularly X-2p states. Figure [5] shows 
the atomic 2s and 2p levels for each of X from B to F 
together with the spin-unpolarized 3d and 4s levels of Fe. 
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FIG. 7. (Color online) Partial density of states of NdFenTlB 
(red solid line) and NdFeiiTii?x (black broken line). The B- 
s, Fe{8j)-s, -p, -d components are shown. The origin of the 
energy is set at the Fermi level. The split-off states can be seen 
around —8eV in the B-s component. Note that the weight 
around +6eV in the Fe(8j) components is not the antibonding 
states coupled with B-s states, but the antibonding states 
coupled with B-p states. 


Although 2s levels are deeper than the Fe-Sd level for all 
X, 2p levels of B and C are definitely higher than the Fe- 
3d level and the N-2p level is nearly degenerate with the 
Fe-3d level as pointed out by Kanamori.— In order to an¬ 
alyze the nature of the X-2p states, most of the available 
theoretical papers take account of only the hybridization 
between X-2p states and Fe-3d states. However, it is 
important to take account of the Fe-4s, p states which 
form broad sp bands overlapping with the Fe-3d bands 
in NdFeiiTiA like in ordinary transition metals^^Z In the 
following discussion, we use the local coordinate as shown 
in Fig. [51 The z axis is parallel to the crystal c axis. 

The B-s and Fe(8j)-s, -p, -d components of the partial 
DOS for NdFeiiTiB are shown in Fig. [71 (For other X 
elements, see Figs. S2 and S3 of the Supplemental Ma¬ 
terials^.) We take a projection of the eigenstates within 
an atomic sphere with the radius Tc- (See Fig. SI of the 
Supplemental Material.) A split-off band at around —8 
eV is present in the B-s DOS and a significant weight ex¬ 
ists also near the bottom of the band above —7 eV. Such 
a feature is correlated with the feature existing in the 
s-component DOS of the Fe(8j)s which are the atoms 
closest to X. [Note that all partial DOSs associated with 
Fe(8j) are for four Fe(8j) atoms. The local x axis is in 
the direction from X to each Fe(8j) atom.] There is less 
significant weight of the p and d states of Fe(8j) in the 
energy range of the B-2s band. Therefore, the states in 
the split-off B-2s band and near the bottom of the con¬ 
tinuous band are dominantly formed by the hybridization 
between B-2s and Fe-4s states. It is also important to 
note that there is no clear evidence of the antibonding 
states between B-2s and Fe-4s states pushed up above 


the Fermi level. This implies that one state is added as 
the occupied state for each spin state by introducing B 
to NdFeiiTi. The split-off X-2s band becomes deeper 
for higher valence elements, C to F, and the 2s band in¬ 
troduces one additional occupied state per spin for all X. 
This is an important result in the later discussion. 

The character of the states in NdFenTiA associated 
with the X-2p states needs more careful analysis. The 
overall features of the partial DOS for X-2p states across 
the typical elements X are shown in Figs. S2 and S3 
in the Supplemental Material^^. For X-2px and 2py 
states, both a- and 7r-type hybridizations are possible 
with Fe(8_)) states. We begin the discussion with the a- 
type hybridization because it is stronger. For the pair 
of X and Fe(8j)o located along the x direction, the 2px 
state of X and s, Px, d^z'^-r'^ and states of the 

Fe( 8 j)o are involved (Fig. [ 8 ]). The sharp structures seen 
in the X-2px DOS are reflected in the DOSs of these 
states of Fe( 8 j) more or less depending on the strength 
of the hybridization. However, we see that the structure 
at the band bottom can also be seen in the DOSs of Py 
and dxy of Fe( 8 j) which do not form cr-type hybridiza¬ 
tion with X-2px. Such a band bottom structure cannot 
be understood in terms of 7r-type hybridization because 
it is even stronger for dxy than for dx'^-y^. The domi¬ 
nant mechanism of producing the structure is due to the 
following ( 7 -type hybridization path; i.e., 2py at X hy¬ 
bridizes strongly with s of Fe( 8 j)i and Fe( 8 j )3 of Fig. [ 6 ] 
and then these s states hybridize with py and dxy states 
of Fe(8j)o. On the other hand, the structures in the X- 
2pz component is of 7r-type origin and the s states at 
Fe(8j) do not take part in the hybridization because four 
Fe(8j) atoms and the X atom form a plane perpendicu¬ 
lar to the crystal c axis which is taken as the z-direction. 
Fe{Sj)o-dzx states and also the Fe(8j)o-Pz state form the 
TT-type hybridization with the X-2pz state. 

From Figs. |5| and ini one may think that the X-2p com¬ 
ponents near the bottom of the host band including the 
split-off bands may be characterized as the X-2p dom¬ 
inating states. However, as pointed out above, for B 
and C at least, the 2p levels are located above the Fe-3d 
levels. Therefore, the structures near the band bottom 
correspond to the bonding states among X-2p, Fe-3d, 
and (for cr) Fe-4s,p states and the weight of B-2p com¬ 
ponent near the band bottom is small as can be seen in 
Figs. |5| and ini In this energy range, both the Fe(8j)-3d 
and 4s components for the cr type and only Fe(8j)-3d 
components for the tt type have larger weight. As we 
move from B to F, the partial DOS of the X-2p compo¬ 
nent shifts to a lower energy and the weight increases. 
The difference between cr type and tt type in the energy 
region of the bonding states comes from the difference 
in the hybridization strength. The pda hybridization is 
stronger than the pdn hybridization by a/S in the canon¬ 
ical band picture^^ and moreover Fe(8j)-4s states do not 
contribute to the rr-type hybridization. 

The behavior of the antibonding states is quite differ¬ 
ent between the a bond and the tt bond. For the a bond. 
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FIG. 8. (Color online) The partial DOS relating to the a bond between X-p and Fe{8j) states [X-p^, Py, Fe{8j)-s, Px, d^^ 2 _^ 2 , 
d^ 2 _y 2 ]. Moreover, Fe(8ji)-pj, and d^y are also shown. The red solid lines are for NdFenTiX and the black broken lines are for 
NdFeiiTlFlx. The origin of the energy is set at the Fermi level. 


significant weight is present a bit higher than 6 eV from 
the Fermi level for B and this weight comes down to a 
lower energy for C and N. The structure becomes weak 
for O and nearly disappears for F. This is due to the 
deeper 2p levels and shrinking of the 2p wave functions 
for O and F. Because of the involvement of Fe(8j)-4s 
states in the cr bond, continuous spectra exist between 


the sharp structures of bonding and antibonding states. 
For B and C, however, the weight of the continuous spec¬ 
tra below the Fermi level is small and only the bonding 
states are occupied. Moving to N and O, we see that 
the weight of the continuous spectra particularly in the 
majority-spin band increases slightly for the X-2px^ 2py 
states even below the Fermi level. Therefore, not only 




































































FIG. 9. (Color online) The partial DOS relating to the tt bond between X-p and Fe(8j) states [X-pz, Fe{8j)-pz, dzx]- Nd-p^ 
and d^^ 2 _z .2 are also shown. The red solid lines are for NdFeiiTiX and the black broken lines are for NdFeiiTii?x. The origin 
of the energy is set at the Fermi level. 


the bonding states but also some part of the antibonding 
states are occupied. The situation is different for the tt 
bond. As Fe(8j)-4s states are not involved in this case, 
the antibonding states are just above the Fermi level al¬ 
ready for B, about to be occupied for C, and occupied in 
the majority-spin band for N. Up to N, the tt antibonding 
state is above the Fe d band and sharp. However, for O 
and F, the antibonding state moves into the d band and 
becomes broad and small. The antibonding states start 
to be occupied also in the minority-spin state for O and 
even more occupied for F. This behavior of the antibond¬ 
ing TT state can be clearly seen in X-p^. and Fe(8j)-dzx 
in Fig. ini Comparing these two DOSs, we find that the 
main component of the tt antibonding state is Fe(8j)- 
dzx- Traces of the antibonding tt state can also be seen 
for Nd-pz, d 322 _j .2 states as shown in Fig. [9]. The differ¬ 
ence between C and N in the filling of the antibonding 
states in the majority-spin band was pointed out before 


for i?Fei 2 A— and also for i? 2 Fei 7 A 3 .— We note that 
in both RFei 2 X and i? 2 Fei 7 A 3 , X and its nearest neigh¬ 
bor Fe atoms form a plane perpendicular to the direction 
from X to R and that the s states at the Fe sites do 
not hybridize with the X-2pz state. Therefore, the fill¬ 
ing of the Ptt antibonding states occurs for A=N in both 
systems. 

We make brief comments on the role of Nd in the X- 
2p related states. We find some clear structures related 
to the bonding and antibonding X-2p states but their 
weight is quite small. Based on this observation, we con¬ 
clude that Nd does not play an important role in the 
bonding with the interstitial element X and that the oc¬ 
cupied antibonding state just below the Fermi level for 
N is ascribed to the hybridization between X-2pz and 
Fe{8j)-3dzx orbitals. Note, however, that the hybridiza¬ 
tion between X-2pz and Nd-6pz: ^d^z'^-r'^ orbitals plays 
an important role in the crystal field parameter A 2 as 
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described later. 


B. Optimized Structures 


shrink in space. These features weaken further the X- 
Fe(8j) covalent bond and lead to further increase of the 
ab area. 


C. Magnetic Moment 



FIG. 10. (Color online) ab area of NdFenTiX (X=B, C, N, 
O, F) and NdFeiiTi. 


As TableUshows, the volume expands with the intersti¬ 
tial X, which leads to appreciable magnetovolume effects 
as shown later. We also observe the following interesting 
aspect in Table H] which reflects the bonding character 
related to the X-2p states discussed above. Figure [TUI 
shows the X dependence of the ab area. The broken hor¬ 
izontal line in this figure is the ab area for NdFenTi. 
For X=B, although the volume expands by about 5 %, 
the ab area decreases very slightly. The ab area shrinks 
farther for X=C. Then as X moves from C to F, the 
ab area starts to increase. This trend is in clear con¬ 
trast to the trend of the atomic radius which decreases 
monotonically from 0.88 A for B to 0.58 A for FA^ The 
behavior of the ab area of Fig. [TOl reflects the strength of 
the covalent bond between X-2p orbitals and Fe(8j)-3d, 
4s orbitals. As was explained in the preceding section, 
only their bonding states are occupied for X=B and C. 
Therefore, the covalent bond between B or C and Fe(8j) 
tries to keep the B, C-Fe(8j) distance short. Moving from 
B to C, we notice that the bonding states become deeper 
and that the weight at C increases. These features con¬ 
tribute to a stronger covalent bond between C and Fe(8j) 
because the antibonding states are still mostly unoccu¬ 
pied. The situations for B and C qualitatively explain 
the variation of the ab area in Fig. [TUI On the other 
hand, for Ar=N, the antibonding 2pz state becomes oc¬ 
cupied in the majority-spin band. At the same time, the 
tail of the cr-type antibonding states are also partly occu¬ 
pied. Therefore the covalent bond between N and Fe(8j) 
is weakened and the ab area becomes larger. As we move 
to O and F, the antibonding states will be further oc¬ 
cupied, 2p levels become much deeper, and 2p orbitals 



FIG. 11. (Golor online) The magnetic moment per formula 
unit for NdFeiiTiX (X=B, G, N, O, F) and NdFenTi. 



FIG. 12. (Golor online) The difference of the spin-resolved 
number of electrons between NdFenTiX and NdFenTiAx- 
Those of the majority-spin states and minority-spin states are 
shown as the filled red circles and blue squares, respectively. 

We show by solid circles in Fig. |TT] the magnetic mo¬ 
ment per formula unit for NdFenTiAl for the optimized 
structures. The broken line near the bottom of the figure 
shows the magnetic moment of NdFenTi with its opti¬ 
mized structure and the open circles near the middle of 
the figure show the magnetic moment of NdFenTiiJjjs:. 
The change in the magnetic moment from the bottom 
broken line to the open circles corresponds to the mag- 
netovolume effect, which is associated with the expanded 
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FIG. 13. (Color online) A schematic picture for explaining the behavior of AA^-|-(X) for A=B and C. 
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FIG. 14. (Color online) The number of local majority-spin (left) and minority-spin electrons (middle) and the spin moment 
(right) at Fe sites 8j, 8f, and 8i of NdFeiiTiX and NdFeiiTiFix. The integration region is given by a sphere with a radius 
giving the same volume of each Bader region. 


volume of NdFenTiiSx compared with that of NdFenTi 
by 4 to 7 % depending on X. However, among the open 
circles, the variation in the magnetic moment does not re¬ 
flect the variation in the unit cell size. On the other hand, 
the change in the magnetic moment from the open cir¬ 
cles to the solid circles corresponds to the chemical effect, 
which is the main subject in this section. Such a mag¬ 
netovolume effect and a chemical effect in materials such 
as Fe 4 N^, and 7 ^ 2 Fe 17 X 3 ^ 1^125 

systems were discussed before. 

In order to explain the trend in the chemical effect on 
the magnetic moment change, we show in Fig. [T^ the 


change in the number of occupied states per unit cell 
for the spin state cr, ANcr{X), caused by adding X into 
NdFeiiTii?x- In this figure and also in the following 
argument, we use the following convention. As already 
pointed out, the 2s states of X form split-off bands be¬ 
low the host valence band and introduce two additional 
occupied states to the host band. Therefore, we assume 
the valence of X, Zx, as 1, 2, 3, 4, and 5 for B, C, N, 
O, and F, respectively, by focusing on only 2p states of 
X. Then ANcr{X) should satisfy the charge neutrality 
condition given by 

AN^iX) + AN^iX) = Zx, 


( 4 ) 
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On the other hand, the chemical effect part of the mag¬ 
netic moment, Am, in Fig. [TT]is given by 

Am{X) = AN-f{X) - AN^{X). (5) 

There are a couple of important aspects in Fig. [T^l 
First, AN^{X) is quite small for X=B and C. This can 
be understood in the following way. To simplify the sit¬ 
uation, we tentatively neglect the presence of wide con¬ 
tinuous 4s, 4p bands above the Fe-3d bands and assume 
that the Fermi level in the majority-spin state is located 
in an energy region where there are no states. Such sim¬ 
plification is a good approximation as long as the anti¬ 
bonding states do not make any significant contributions 
to the occupied states. Moreover, for the states, the 
contribution from the wide 4s, 4p bands is quite small. 
Figure IlSa) shows schematically the DOS of the host 
material. The Fermi level of the majority-spin state, Ep, 
is located just above the host band. As shown in Fig. 0 
B and C have their 2p levels located above the Fe-Sd 
band. With introducing the hybridization between the 
host states and the X-2p state, a bonding state appears 
below the host band and the antibonding state above it 
without being occupied. The bonding and antibonding 
states are formed from the original host states and the 
X-2p state. The total number of the host states used to 
form the bonding and antibonding states is just one as 
long as the weight of the X-2p state is one. Therefore, 
for the majority-spin state where the host band is fully 
occupied, while one additional state is introduced as the 
bonding state, the same number of states is subtracted 
from the host bands keeping the number of occupied 
states unchanged. This situation is shown in Figs. fT^ dl 
and fT^ el. Figure IT^dl shows the difference in the DOS 
between Fig. [TSjb) and Fig. fT^ al. The bonding and 
antibonding states are additionally introduced but some 
number of states whose total weight is one is subtracted 
from the host band. Figure me) shows the integrated 
DOS obtained by integrating the DOS of Fig. fT^dl over 
energy. This quantity jumps up to unity when the up¬ 
per limit of the integral crosses the bonding level, then 
starts to decrease as the upper limit of the integral moves 
into the host band region and reaches zero above the 
host band. It jumps up to unity again when the up¬ 
per limit of the integral crosses the antibonding state. 
Therefore, as long as the antibonding state is unoccu¬ 
pied, the change in the number of occupied states per 
unit cell for the majority-spin state AA^'|-(A) is zero. At 
the same time. Fig. me) also shows that in the minority- 
spin band where the Fermi level is located within the host 
band, ANi{X) > 0. With the above simplified model, we 
obtain that AN^{X) = 0 for W = B, C and ANi(X)=l 
for X=B and 2 for X=C to satisfy Eq. (4). Accordingly, 
Eq. (5) gives as the chemical effect part of the magnetic 
moment Am{X) = —1 for X—B and —2 for X=C. These 
results are qualitatively consistent with the behavior seen 
in Eigs. ITT] and \n\ for X=B and C. The deviation of the 
results in these figures from those of the simplified model 
comes from the partial filling of the antibonding states 


through the wide continuous 4s and 4p bands above the 
host band. 


Second, the sudden increase in AN^{X) from C to N 
is due to the filling of the antibonding state for the 
majority-spin state (Eig. [5]) and leads to the sudden in¬ 
crease in Am{X) in Eig.[TT] Eurther increase of AN^{X) 
from X=N to E is due to the filling of p^ antibonding 
states. Third, compared with the behavior of AN^{X), 
the variation of ANi{X) is small. This is partly due to 
the fact that the Eermi level of the minority-spin state is 
located at the deep valley in the bulk density of states. 
To stabilize the band energy, the Eermi level is pinned 
in the valley region. Another important reason for the 
small variation of AN{Ep, X) is that the Eermi level of 
the minority-spin state is deep inside the d band. This 
makes it difficult for the antibonding state to be occu¬ 
pied. Therefore, the minority-spin antibonding state 
is filled only for A=F though its weight is small due to 
the reduced pdir hybridization mentioned before. 


We now move to the discussion on the variation of 
the magnetic moment of Ee among 8j, 8/, and 8i sites, 
which is shown in Eig. [TH Because of the presence of 
Ti at one of the 8i sites, the electronic structure is not 
the same among Fe atoms even within a given category 
of sites. Therefore, the average value is shown for each 
of 8j, 8/, and 8* sites. Kanamori explained the behavior 
seen in Fig. [T3] using the concept of cobaltization^ It 
says that by introducing B or C, the hybridization re¬ 
pulsion between X-2p and Fe-Sd states pushes down the 
minority-spin 3d states of Fe(8j) more strongly than its 
majority-spin states leading to significant increase in the 
minority-spin electrons. Slight decrease in the majority- 
spin electrons is a result of the balance between the hy¬ 
bridization repulsion and the local charge neutrality re¬ 
quirement in metals, with the Fe(8j) being slightly more 
electron populated. The resulting electronic structure of 
Fe(8j) is similar to that of Co and the magnetic moment 
of Fe(8j) is reduced by about 0.4 pe- On the other hand, 
the change in the electronic structure at Fe(8/) is just op¬ 
posite to that at Fe(8j). Kanamori also pointed out that 
the behavior at Fe(8/) can be explained in parallel with 
the change at the Fe sites next to Co in the Fe-Co alloy. 
The above picture of cobaltization caused by the inter¬ 
stitial B is schematically explained by Figs. 1 and 2 in 
one of his papersi^ The concept of cobaltization seems 
to be reasonable for X=B and C. However, the net mag¬ 
netic moment is decreased by the chemical effect, and 
the enhancement by X is due to the magnetovolume ef¬ 
fect. For X=N, O, and F, the atomic 2p level is very 
close to or below the Fe-3d level (Fig. [S]) and the picture 
of cobaltization is not applicable. The positive chemical 
effect for these elements is due mainly to coming down 
of the Ptt antibonding state below the Fermi level in the 
majority-spin state. 
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FIG. 15. (Color online) (r^)j42 of NdFenTiX. Contributions 
from the static Coulomb and the exchange-correlation poten¬ 
tial are shown by black circles and green circles, respectively. 
Red circles show the sum of the two contributions. 


FIG. 17. (Color online) Static Coulomb contribution to 
{r^)A 2 brought by the off-site charge density (outside the Nd 
site). 



FIG. 16. (Color online) p-p, d-d, and s-d components in 
(r^)j42. The cutoff radii are determined from Bader analy¬ 
sis. 


the contribution from the static Coulomb potential. 

The static Coulomb potential is caused by the charge 
density. We separate it into the on-site charge (at the 
Nd site) and off-site charge. Since the Nd-ion charge is 
spherically symmetric, the on-site contribution to 
is a consequence of the aspherical (e.g., p, d) compo¬ 
nent of the electron density at the Nd-site. To see this 
in more detail, we decompose the electron density within 
the Bader sphere at the Nd-site into the components hav¬ 
ing as follows^. 


V-a{r) = 
p{R) = 


f 

Ir<r„ k ~ -R| 


( 6 ) 

( 7 ) 


-E E E fummiR) 

i I m— — l 

- Y. ne,imi>m>{R)Yi^{RyYi,^,{R) ( 8 ) 

Iml'm' 


where. 


D. Crystal field parameter 


n,,lml'myR) = Yflrn(Rrfi'n.'(R) ( 9 ) 


Finally, we discuss the origin of the dopant dependence 
in {r^A^. In order to see whether the dopant depen¬ 
dence can be explained from the structure change, we 
plot {r^A^ of NdFeiiTii?x in Fig. [3] by open circles. 
The difference among X=B, C, N, O, and F is within 
about 220 K and the structural effect is not sufficient to 
explain the dopant dependence. 

The potential used in [Eq. ([3])] consists of the 

static Coulomb (Flartree) and the exchange-correlation 
potentials. Figure [15] shows the dopant dependence of 
each contribution. We can see that the static Coulomb 
potential dominantly affects the dopant dependence of 
(r^)A 2 . In the succeeding analysis, we concentrate on 


//m(77) is the radial part of the Im component in the ith 
eigenstate. The cutoff radius Rc of the integration in 
Eq. ([3|) is determined from the Bader region mentioned 
above (=rc). [We note that Rc is conceptually different 
from Tc, although the same value is taken in the actual 
calculations. The former is the cutoff radius of the charge 
density to evaluate the static Coulomb potential, whereas 
the latter is the cutoff radius for the potential in Eq.®.] 
l/|r — i2| can be expanded by the real spherical harmon¬ 
ics: 


1 


E 


47r 

+ 


Zr{r)Zy{R) 


( 10 ) 
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FIG. 18. (Color online) pz and d^^ 2_^2 components (rie.p^ = ne.ioio and ne,d 2 _ 2 = ^£, 2020 ) of the electron density difference 
between NdFenTiX and NdFenTiFix are shown on the left and right panels, respectively. The horizontal axis is the distance 
from Nd. The positive value means that is increased by introducing X. The black arrows indicate the Bader radius. 



FIG. 19. (Color online) Electron density difference between NdFenTiX and NdFeiiTiEx- The density increases (decreases) 
in the red (blue) region. The figures show the density difference map on the ac plane around the Nd atom. The Nd atom is at 
the center of each figure. The vertical and horizontal ranges are 0.2 times c and a axis, respectively. The displayed region is 
indicated by a black dashed box on the right panel. The outermost red island corresponds to the peak at 0.25A in Fig. 1181 


Here, r< = min(r, i?) and r> = max(r, i?). In¬ 
serting Eqs. (0), ®, and m' into Eq.(j6l), and per¬ 
forming angular integration of ZY^{R)Yli^ {R 
three components, 

and rie ,0020 + i^e, 2000 , yield nonzero contribution to A^. 
The dopant dependencies of the three components, p-p 

{— — 1 d-d (= — 2 ^e,2m2m)7 and S- 

d (= nepo 2 o + R'e, 20 oo)j ^re shown in Eig. [121 Among 
the three, the p-p electron density gives a peak structure 
at X=N, which is similar to the behavior of the total 
In contrast, the contribution from d-d density 
decreases monotonically from X =B to F. The s-d cross¬ 
term contribution is small compared with the p-p and d-d 
contributions. We also plot the remaining component de¬ 
fined as the static Coulomb minus the contribution from 
the electron density within the Bader atomic sphere in 
Fig. [T71 The contribution increases monotonically as the 
dopant changes from B to F. 

In order to understand the dopant dependence of the 
p-p and d-d components, we decompose them into each m 
component, and compare the results for NdFenTiAl and 
for NdFeiiTiEjf. We found that most of the components 
are not changed between NdFenTiX and NdFeiiTiEjf. 
However, the p^ and d^z 2_^2 components are significantly 
changed as shown in Fig.jTHl where the ne,imi'm'{R) dif¬ 
ferences for Pz and d^z^-r^ are plotted as a function of 


the distance R from Nd. The changes are prominent 
in the vicinity of the Nd atom within a 0.5 A radius. 
In Fig. [T21 the density difference map on the ac plane 
including the Nd atom is shown. Both pz and d^j 2_^2 
orbitals have an anisotropic shape expanding along the 
c axis. Correspondingly, an increase in the density dif¬ 
ference is observed above and below the Nd site along 
the c axis. This accumulated electron charge suppresses 
the aspherical 4/-electron density from above and be¬ 
low the Nd site, which induces the uniaxial anisotropy. 
For Pz density, the difference increases from X—^ to N 
and decreases from N to F, and for d^z^-r^, it decreases 
monotonically from X—^ to F. These density differences 
near the Nd site explain the dopant dependence of the p- 
and d-components. The influence of the bonding charge 
between X and R has also been discussed in Ref. HO- 
There still remains the contribution in the static 
Coulomb potential other than the contribution from the 
electron density within the Bader atomic sphere. This 
can be interpreted as the off-site contribution, namely 
the contribution from the charges outside the Nd-atomic 
sphere with Re- The increasing behavior of this con¬ 
tribution shown in Fig. [T7] reflects the electronegativity 
of the interstitial dopants. The dopant elements B, C, 
N, O, and F have the electronegativity 2.04, 2.55, 3.04, 
3.44, 3.98, respectivelyj^ F tends to attract the negative 
charge more than B. The attracted negative charge at 
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the X site gives the Coulomb repulsive interaction at the 
Nd site along the c axis. This contribution for 
has a positive value. Thus, F gives a less negative value 
than B. 

These three contributions, p-p, d-d electron charges at 
Nd-site, and the charge at the interstitial dopant, explain 
the nonmonotonic dopant dependence of 

V. CONCLUSION 

We have studied the magnetic properties of 
NdFeiiTiX with X=B, C, N, O, and F with first- 
principles electronic structure calculations to find a 
better interstitial dopant for a permanent magnet. 
According to our study, the total magnetic moment is 
increased by doping B, C, N, O, and F at the interstitial 
site. The magnetic moment increase can be separated 
into the magnetovolume effect and the chemical effect, 
the latter of which depends strongly on the dopant. The 
chemical effect is negative for B and C and becomes 
suddenly positive for N, O, and F. By studying carefully 
the basic electronic structures of NdFenTiA from the 
viewpoint of the hybridization between X and Fe(8j), 
we gave an explanation to the mechanism of the negative 
chemical effect for B and C. The sudden increase in the 
chemical effect part of the magnetic moment from C to 
N is due to the coming down below the Fermi level of 
the majority-spin p^^. antibonding state formed between 


X and Fe(8j). 

The results for {r‘^)A^ suggest that the strongest en¬ 
hancement of uniaxial anisotropy is achieved by A=N. 
We analyzed the dopant dependence of from the 

decomposed electron density ne,imi'm'{R)- The analysis 
suggests that the bonding charge between Nd and A, and 
the charge at the X site can explain the increasing trend 
from B to N and decreasing trend from N to F of the 
magnetocrystalline anisotropy. 

The present study suggests that the interstitial nitro- 
genation is the most appropriate among the typical ele¬ 
ments B, C, N, O, and F doping in terms of the magne¬ 
tization and magnetocrystalline anisotropy. 
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VI. SUPPLEMENTAL MATERIAL 

In this Supplemental Material, we show the following 
figures. 

1. Fig. ISlI The radius of the atomic sphere whose vol¬ 
ume is that of the Bader region^ii^ for each atom 
in the unit cell of NdFenTiA. The same atomic 
sphere radius is used for NdFenTiFlx. 

2. Fig. [52] The partial DOS of X in NdFenTiX. 

3. Fig. jSSj The zoomed figures of Fig. |S2l 

4. Fig.jSH The partial DOS of Fe(8j) in NdFenTiX. 

5. Fig.jSSj The partial DOS of Nd in NdFenTiX. 

6 . Fig.|S6| The number of electrons for each spin state 
in the atomic sphere for each X in NdFenTiX. 

The partial DOSs in Figs. IS2[ IS3[ IS4I and |S5] and the 
number of electrons in Fig. IS6I are given for the atomic 
sphere with the radius given in Fig. jSTj X-s and Px +Py 
components in Figs. |S2] and |S2] are related with the cr 
bond between X and Fe(8j), while X-p^ component is 
related with the tt bond. In Fig. |S6l the X-p^ compo¬ 
nent shows a jump from X=C to N. This is due to the 
occupation of the anti-bonding states at X=N. 


atomic radius [ A ] 
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FIG. SI. (Color online) The sphere radius used for the inte¬ 
gration region of the partial DOS for each X. 
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FIG. S2. (Color online) Partial density of states at X of NdFenTiX. The s components (upper panel), the sumation of Px 
and Py components (middle panel), and pz components (lower panel) are shown. 



FIG. S3. (Color online) Zoomed figures for partial density of states at X of NdFenTiX corresponding to Fig. IS2I 
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FIG. S4. (Color online) Partial density of states at Fe(8y) of NdFenTiX. The s, Px, Py, Pz, d^y, dyz, d^^2_x2, 
components are shown. The red solid lines for NdFenTiX and the black broken lines for NdFenTiFlx. 
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FIG. S5. (Color online) Partial density of states at Nd of NdFenTiX. The s, p^, Py, Pz, dxy, dyz, d^^ 2 _x 2 , 
components are shown. The red solid lines for NdFenTiX and the black broken lines for NdFenTiiJx. 
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B C N O F 


FIG. S6. (Color online) The number of majority spin and 
that of minority spin electrons within the sphere around X. 
The components are seperated into each orbitals, s, Px +Py, 
and Pz ■ 






